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Abstract 

Species richness on oceanic islands has been related to a series of ecological factors including island size and isolation (i.e. 
the Equilibrium Model of Island Biogeography, EMIB), habitat diversity, climate (i.e., temperature and precipitation) and 
more recently island ontogeny (i.e. the General Dynamic Model of oceanic island biogeography, GDM). Here we evaluate the 
relationship of these factors with the diversity of bryophytes in the Macaronesian region (Azores, Madeira, Canary Islands 
and Cape Verde). The predictive power of EMIB, habitat diversity, climate and the GDM on total bryophyte richness, as well 
as moss and liverwort richness (the two dominant bryophyte groups), was evaluated through ordinary least squares 
regressions. After choosing the best subset of variables using inference statistics, we used partial regression analyses to 
identify the independent and shared effects of each model. The variables included within each model were similar for 
mosses and liverworts, with orographic mist layer being one of the most important predictors of richness. Models 
combining climate with either the GDM or habitat diversity explained most of richness variation (up to 91%). There was a 
high portion of shared variance between all pairwise combinations of factors in mosses, while in liverworts around half of 
the variability in species richness was accounted for exclusively by climate. Our results suggest that the effects of climate 
and habitat are strong and prevalent in this region, while geographical factors have limited influence on Macaronesian 
bryophyte diversity. Although climate is of great importance for liverwort richness, in mosses its effect is similar to or, at 
least, indiscernible from the effect of habitat diversity and, strikingly, the effect of island ontogeny. These results indicate 
that for highly vagile taxa on oceanic islands, the dispersal process may be less important for successful colonization than 
the availability of suitable ecological conditions during the establishment phase. 
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Introduction 

The Equilibrium Model of Island Biogeography (EMIB) states 
that, other things being equal, area and geographic isolation are 
the two main factors determining extinction and immigration 
rates, which in turn regulate the level of species richness that is 
reached at a dynamic equilibrium [1], [2]. Although many 
hypotheses have been proposed to explain the role of area in 
species richness patterns [3], in its original formulation the EMIB 
postulated the effect of area per se referring specifically to 
demographic processes (i.e. smaller areas support smaller popu- 
lations that are hence more prone to species extinctions). Despite 



its importance in the development of ecology and biogeography, 
the EMIB has been criticized for the lack of ability of its simple 
mechanisms to account for variations in species richness (e.g. [4]). 
In fact, models based on additional factors have been suggested to 
also account for island diversity, including energy [5] , [6] , habitat 
diversity [7] or island ontogeny in the particular case of oceanic 
archipelagos [8]. 

In essence, the models considering energy relate the amount of 
available resources with the possibility of maintaining higher 
population sizes and therefore more species [6]. The variety of 
resource types (e.g. habitat diversity) would also promote the 



PLOS ONE | www.plosone.org 



1 



July 2014 | Volume 9 | Issue 7 | e101786 



Determinants of Bryophyte Richness in the Macaronesia 



coexistence of more species by diminishing interspecific compe- 
tition and increasing sympatric speciation through ecological space 
partitioning (see [9]). The ontogenetic evolution of the island itself 
may affect its carrying capacity and hence species richness, 
because the variations in area and structural complexity occurring 
during the island's life cycle influence both typical immigration- 
extinction dynamics and diversification by in situ speciation (i.e. the 
General Dynamic Model of oceanic island biogeography or GDM) 
[8], [10]. The effects of energy, habitat diversity and island 
ontogeny on species richness have been typically examined using 
surrogates such as actual evapotranspiration or other climatic 
factors (e.g. [1 1]), topographic variables or habitat classifications 
(e.g. [12]), and the maximum geological age dated for islands (e.g. 
[13], [14]), respectively. Numerous studies have evaluated some of 
these models for a wide variety of taxa and archipelagos, either 
confirming or rejecting their predictions (e.g. [15] — [17]). 

Although all these factors are known to affect island species 
richness, few attempts have been made to assess their comparative 
importance within a single evaluation (but see [11], [18], [19]). 
This may be due to the fact that most predictors are often 
correlated and therefore it is difficult to separate their true 
influence on species richness through common statistical tech- 
niques [20] (see also [21]). In addition, generalizations about the 
importance of the processes underlying these predictors depend on 
the idiosyncratic characteristics of both islands (e.g. the range of 
variation in area, isolation or elevation) and taxa (e.g. dispersal 
ability or life cycle). For example, the influence of isolation on 
immigration depends on the dispersal ability of the taxon, which in 
turn limits the probability of in situ speciation [14], [22], [23]. 
Similarly, the influence of environmental heterogeneity (habitat or 
climatic diversity) on the successful establishment of species varies 
according to their physiological and ecological tolerances (i.e. 
niche breadth), eventually determining the shape of richness- 
environment relationships (see [24]). 

Bryophytes - which encompass hornworts, liverworts and 
mosses - are unique among land plants because: (i) the 
gametophyte is the dominant phase of the life cycle comprising 
the leafy or thalloid plants; and (ii) the sporophyte, which consists 
mainly of a short-lived small "capsule", is always attached to and 
dependent on the gametophyte. Both singularities make the two 
generations of the life cycle to contribute significantiy to the 
dispersal and establishment processes [25] . In addition, contrary to 
seed plants and ferns, they lack complex vascular tissues and 
developed a poikilohydric strategy that allows them to absorb 
water over their whole surface by capillarity, being able to remain 
metabolically inactive when dry conditions exist. Furthermore, 
bryophytes are characterized by extremely low levels of endemism 
in oceanic floras (see [26] for review], which is thought to be a 
consequence of the high dispersal ability of the group [27]. 

Despite these interesting features, bryophytes have received 
relatively little attention in island biogeography studies compared 
to other plant groups (but see [22], [28]-[30]). Most of these works 
include only one archipelago (but see [17], [22]) or do not consider 
all the above-mentioned factors, and in particular climate (but see 
[29]). Also, the effect of climate on large-scale species richness 
gradients has been occasionally analyzed in spore-dispersed plants 
[31], being mostiy studied indirectly through its correlation with 
latitude and altitude (e.g. [32]-[34]). In the present study we 
examine the role of geographical, temporal and environmental 
factors on the between-island variation of bryophyte species 
richness in the Macaronesian Region (i.e. Azores, Madeira, 
Canary Islands and Cape Verde). Specifically, we evaluate four 
non-exclusive hypotheses under the following premises: 



H t . The Equilibrium Model of Island Biogeography (EMIB) 
should not significantly account for the variation in species 
richness of bryophytes, or its effect should be negligible. We expect 
that geographic isolation will not have a significant effect on 
immigration rates since bryophytes have the potential to disperse 
long distances by spores [35], [36]. In spite of some discrepancies 
[17], [22], [37], the dispersal ability of bryophytes should in turn 
limit the influence of area/w se because the high rescue effect from 
surrounding source populations would minimize species extinc- 
tions. 

H 2 . The General Dynamic Model of oceanic island biogeog- 
raphy (GDM) should not be of high relevance for bryophytes 
because the effect of area in species richness should be minimized 
with increasing dispersal ability and also because former studies 
suggested that time per se appears to have little support in 
predicting species richness in the group [17]. 

H 3 . Habitat diversity (HD) should have a significant effect on 
species richness because bryophyte communities are known to 
show significant degrees of compositional turnover between 
different habitats [29], [38]-[40]. 

H4. Precipitation and temperature (CLIMATE) should have a 
strong effect on species richness, since sexual reproduction and 
photosynthesis in bryophytes are highly dependent on water 
availability, and optimal growth occurs with moderate tempera- 
tures [41]. 

Results 

Univariate regressions between dependent variables (S T0T , S M 
and S L ) and all the considered predictors showed similar results for 
mosses and liverworts (Table 1). From the predictors representing 
the EMIB only area (A) was significantly related to moss species 
richness variation. In the case of GDM, both the linear and 
quadratic functions of time were not statistically significant for any 
of the groups. For the HD hypothesis, however, most variables 
were correlated with species richness of mosses and liverworts, 
being highly significant in the former group. Regarding the 
CLIMATE hypothesis, orographic mist layer (MistL) accounted for 
the highest proportion of data variability in both S M and Si, The 
negative effect of maximum temperature (7a£4x) on liverwort 
species richness was also remarkable. Temperature seasonality 
(Ty), although showing a lower correlation, was statistically 
significant for both mosses and liverworts. 

The subset of variables included in the best model for each 
hypothesis was also similar between mosses and liverworts 
(Table 2). The EMIB never exceeded 22% of explained variance, 
being weakly or even marginally significant, while GDM was 
statistically significant for both groups of bryophytes, particularly 
in mosses. Here note that while time alone is a poor predictor of 
species richness (Table 1), when included in a model with area (^4) 
the two variables account for as much as 64% and 36% of data 
variation in moss and liverwort species richness, respectively 
(Table 2). HD seems to be particularly important for mosses, 
although CLIMATE was the model with the highest explanatory 
capacity for all groups (up to 77%). Given that the above 
mentioned differences between mosses and liverworts cannot be 
discerned when considering all species together (Sj-ot), henceforth 
we will focus on comparing the main findings for both taxonomic 
groups separately. 

Results from partial regressions including the hypotheses that 
seem to better explain species richness (see Table 2) indicated that 
the combined 'GDM+CLIMATE' model explained most of the 
variation in moss and liverwort richness (87.0% and 91.1%, 
respectively), followed by 'HD+CLIMATE' (71.5% and 79.8%) 
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Table 1. Univariate regressions explaining the variation in species richness of all Macaronesian bryophytes (S TOr ), mosses (S M ) and 
liverworts (SJ as a function of the predictors chosen for the Equilibrium Model of Island Biogeography (EMIB), the General Dynamic 
Model (GDM), the Habitat Diversity model (HD) and the Climatic Model (CLIMATE). 







All bryophytes (5 TOT ) 


Mosses (5 M ) 




Liverworts (5 £ ) 




R 2 


F 


R 2 


F 


R 2 


F 


EMIB 


A (+) 


0.132 


2.59 


0.222 


4.85* 


0.015 


0.25 


Dm 


0.002 


0.43 


0.009 


0.16 


0.102 


1.94 


D, 


0.017 


0.30 


0.037 


0.65 


< 0.001 


0.00 


N 


0.079 


1.46 


0.064 


1.16 


0.087 


1.62 


GDM 


T 


0.094 


1.76 


0.047 


0.85 


0.189 


3.97 r 


Tf 


0.191 


1.88 


0.154 


1.46 


0.243 


2.57 


HD 


ELEV (+) 


0.405 


11.58** 


0.482 


15.83*** 


0.201 


4.26* 


sdELEV (+) 


0.505 


17.32*** 


0.606 


26.13*** 


0.242 


5.42* 


SLOPEdiv 


0.045 


0.79 


0.050 


0.89 


0.023 


0.39 


EZ{+) 


0.644 


39.75*** 


0.579 


23.39*** 


0.224 


4.92* 


CLIMATE 


Tmax ( ) 


0.195 


4.1 2 f 


0.082 


1.53 


0.467 


14.87** 


r s (+) 


0.275 


6.45* 


0.234 


5.19* 


0.294 


7.10* 


P MIN 


0.004 


0.07 


0.005 


0.09 


0.103 


1.96 


Ps 


0.006 


0.10 


0.006 


0.10 


0.130 


2.53 


P ANN (+) 


0.033 


0.58 


0.002 


0.03 


0.187 


3.90 1 ' 


MistL (+) 


0.584 


23.90*** 


0.503 


17.20*** 


0.610 


26.62*** 



1 p<0.06, 
*p<0.05, 
**p<0.01, 
***p<0.001. 

Variable codes: A (area), D M (distance to mainland), D, (distance to the nearest island), N (neighbour index); T (oldest geological age); ELEV (maximum elevation), sdELEV 
(standard deviation of elevation), SLOPEdiv (diversity of slopes); EZ (number of ecological zones); T MAX (maximum temperature of warmest month), T s (temperature 
seasonality), P Mm (precipitation of driest quarter), P s (precipitation seasonality), P ANN (annual precipitation), MistL (orographic mist layer). 

The explanatory capacity of each variable {R 2 ) and its statistical significance (F-test) are shown. The sign of the relationship is indicated under parenthesis after the 
predictor variable only when there is a significant effect on the dependent variable. The best fitting function (including significant quadratic functions) is shown in all 
cases, except for GDM for which both linear (T) and quadratic (7T 2 ) functions of time are included as suggested in the literature [8], [10]. 
doi:1 0.1 371 /journal.pone.01 01 786.t001 



and 'HD+GDM' (71.1% and 43.6%) (Fig. 1). However, there were 
contrasting differences in the independent and shared effects of 
these hypotheses between both groups. In mosses, shared effects 
were very high in all pairwise model combinations (ranging 
between 48.9-60.6%), while pure effects were much lower (Fig. 1). 
Hence, when combined with the other hypotheses, CLIMATE 
alone explained around 11.0-19.2% of S M , GDM between 10.5- 
18.9% and HD no more than 3.5%. In the case of liverworts, 
however, around half of the variability in species richness was 
accounted exclusively by CLIMATE, while the independent 
effects of both GDM and HD were weaker (up to 19% and 
0.4% of explained variance, respectively). Note that pairwise 
comparisons between GDM and HD hypotheses showed a 
relatively higher contribution of the former over the latter in both 
taxonomic groups. 

The spatial analyses evidenced that the above models account 
for almost all spatially-structured variation in the data. Model 
residuals were not significantly correlated with latitude in any case 
except for the HD model in liverworts (Spearman r=0.57, 
j& = 0.0 1). In fact, there were only three statistically significant 
autocorrelation values in model residuals (in the case of GDM for 



mosses and of GDM and CLIMATE for liverworts; see Tables 
S3.1, S3. 2 in File S3). Further, archipelago idiosyncrasies seem not 
to have affected model estimates, since SAR models were 
consistent with OLS regressions (File S3) and the statistical 
significance of the parameters for the variables included in OLS 
regressions remained similar. Once the spatial structure was taken 
into account, the predictive ability of SAR models (R z ) increased 
only slightly (compare Table S3. 3 in File S3 with Table 2). 

Discussion 

Although interpreting patterns of species diversity is a recurrent 
issue in island biogeography, the number of studies examining 
different factors altogether is surprisingly low (e.g. [1 1], [18], [42], 
[43]). For instance, Kreft et al. [19] related worldwide patterns of 
vascular plant island species richness to geographic, topographical 
and climatic characteristics, filtering also by the geological origin 
of islands. However, these authors did not account for the different 
intra-archipelago relationships expected when crossing biogeo- 
graphical regions [20] and their meta-analyses did not disentangle 
the combined and independent effect of each factor. Our findings 
suggest that climate and habitat are the most relevant factors in 
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Table 2. Multiple regression results showing the best subset of predictors for each considered model (EMIB, GDM, HD and 
CLIMATE) to explain the between-island variation in species richness of Macaronesian bryophytes. 







F 


P 




AIC C 


Total species richness (5 ror ) 


EMIB (A) 


2.78 


0.126 


0.132 


238.5 


GDM (A, TT 2 ) 


7.93 


0.002 


0.565 


230.1 


HD (sdELEV) 


17.32 


<0.001 


0.565 


227.8 


CLIMATE (T MAX , P MIN , MistL) 


15.67 


<0.001 


0.728 


221.2 


Moss species richness (S M ) 


EMIB (A) 


4.85 


0.042 


0.222 


285.0 


GDM (A, TT 2 ) 


10.54 


<0.001 


0.638 


212.0 


HD (sdELEV) 


26.13 


<0.001 


0.606 


208.8 


CLIMATE (P MN , MistL) 


17.07 


<0.001 


0.662 


208.0 


Liverwort species richness (5J 


EMIB (A, Dm) 


2.85 


0.088 


0.219 


198.8 


GDM (A, TT 2 ) 


3.83 


0.032 


0.363 


197.5 


HD (sdELEV) 


5.42 


0.033 


0.242 


196.1 


CLIMATE (Tmax P Mm MistL) 


19.28 


<0.001 


0.768 


178.3 



The best subset of variables that were chosen using the lowest sample size-corrected Akaike information criterion (AIC C ) is shown in brackets. Adjusted R 2 values and its 
statistical significance according to the F-test are also shown. Model acronyms and variable codes as in Table 1. 
doi:1 0.1 371 /joumal.pone.01 01 786.t002 



The role of climate and habitat in bryophyte diversity 

Our results show that the variation of bryophyte species richness 
between the Macaronesian islands can be relatively well predicted 
by climatic conditions, particularly with those favouring higher 
humidity. In this respect, it is well known that water availability is 
an important component of several key ecophysiological processes 
in bryophytes [25] . One could assume that water availability is of 
major importance during the colonization of oceanic islands by 
bryophytes. Since long-distance dispersal (LDD) is mostiy driven 
by spores (see below), water becomes essential since both sexual 
reproduction and spore germination depends on it. Once the 
populations are effectively established on an island, the availability 
of water at adequate growing temperatures becomes crucial to 
achieve a positive net photosynthetic rate over time, otherwise the 
plants enter dormancy upon drying (reviewed by [44]). Conse- 
quently, the time of survival when an individual colony or shoot 
remains dry is also dependent on the periods with mild 
temperatures. It follows that the role of water availability (and 
secondarily temperature) to maintain populations and ultimately 
avoid local species extinctions could be the main reason to explain 
the relevance of the considered climatic variables in our study. 

Contrary to most vascular plants, bryophytes regulate water 
uptake mainly by capillarity and since they are usually small, their 
survival may be more restricted by the frequency rather than by 
the volume of rainfall [45]. Actually, dew or cloud water 
deposition is often sufficient to remoisten most bryophyte species. 
In agreement with this idea, our results show that orographic mist 
layer (MistL), together with lower values of minimum precipitation 
(Pmin) and maximum temperature (Tmax) correlate with higher 
values of species richness. Such climate is typical of the Azorean 
archipelago and Madeira island but also, to a lesser extent, of La 
Gomera and La Palma in the Canaries [28]. The importance of 
air humidity is clearly evident in (sub-) tropical rainforests [46], 
especially in canopy epiphytes (see [47]), whose species diversity 
may indeed be comparable to some Azorean islands where more 
than 25 species can occur in plots of only 30 cmx30 cm [38]. 



explaining bryophyte richness differences on oceanic islands. 
These suggest that the factors operating during species establish- 
ment may be relatively more important than dispersal during the 
colonization process, at least in the case of the Macaronesian 
region. 




GDM HD HD 

+ + + 

CLIMATE CLIMATE GDM 



Figure 1. Partial regressions showing all pairwise comparisons 
between the models that better fit the species richness 
(Table 2). In all cases it is shown the percentage of variance explained 
exclusively by each model and the shared variance between each pair 
of models. Model acronyms as in Table 1. 
doi:1 0.1 371 /journal.pone.01 01 786.g001 



PLOS ONE | www.plosone.org 



4 



July 2014 | Volume 9 | Issue 7 | e101786 



Determinants of Bryophyte Richness in the Macaronesia 



However, the impact of climate must be broadly analyzed, 
taking into account the landscape conditions where the species 
grow [48] . Our results show high correlations of habitat diversity 
with bryophyte diversity, providing an indirect signal for a certain 
degree of habitat specialization in Macaronesian bryophytes, in 
agreement with evidences that closely related species from this 
group may coexist sympatrically in separate niche spaces [49]. 
This implies that the high dispersal ability of bryophytes does not 
lead necessarily to habitat generalism, as for other taxa with strong 
vagility (e.g. [12]). Strikingly, habitat diversity seems to be as 
important as climate in our case, particularly for the Macaronesian 
moss flora. However in this study, most of the variables 
representing habitat diversity are eminendy topographical and 
hence somehow correlated with island ontogeny (i.e. time) as well 
as with climatic predictors resulting in strong shared effects and 
making it difficult to disentangle their specific influence on species 
richness (see Fig. 1). In fact, the explanatory power of this factor 
remains similar even when we consider the number of ecological 
zones (EQ as a surrogate of habitat diversity, probably because the 
major vegetation formations in Macaronesia are strongly struc- 
tured in altitudinal belts (see File SI); the number of ecological 
zones is highly correlated with both sdELEV (r= 0.94) and ELEV 
(r= 0.92), thus being a surrogate for mesoscale climate gradients as 
well as for habitat diversity per se (see [50]). 

The contrasting patterns found in liverworts and mosses could 
at first be related with their distinct ability to produce sexual and 
asexual diaspores, but no differences in the expression of several 
life-history traits between the two groups were detected in a suite 
of oceanic archipelagos [51]. The apparently higher climatic 
sensitivity of liverworts compared to mosses could then be 
understood by their lower desiccation tolerance [52], especially 
notable in leafy liverworts due to their life-form traits [53]. This is 
mirrored by the high sensitivity of the group to human-induced 
disturbances [39]. Perhaps due to this, mosses can be found in a 
comparatively wider range of landscapes, including grasslands and 
other man-made habitats [25], while liverworts seem to be more 
dependent on sheltered habitats like forests as compared with 
other open landscapes. In fact, liverwort richness at the Azorean 
native forests is higher than that of mosses above 600 m a.s.l. [38]. 
Further studies are necessary to confirm such higher habitat 
specificity across the latitudinal gradient provided by the 
Macaronesian archipelagos. 

Effects of island isolation, area and time 

Long-distance dispersal (LDD) is a rare and stochastic event 
[54], although its prevalence over long time periods may be 
common [4]. In the case of bryophytes, LDD may occur only 
occasionally because spore production seems to be highly 
constrained due to unsuccessful sexual reproduction and even 
when occurring, spore release typically falls within the first tens of 
meters [55] (but see [56]). Hence, asexual propagation is often the 
most frequent way of dispersion [35]. This argument is also 
supported by an increasing number of molecular evidences 
showing that population connectivity at local or even landscape 
scales is hindered by dispersal limitation (e.g. [57], [58]). In line 
with this rationale, significant shifts in life-history traits towards 
decreased sporophyte production and increased production of 
specialized asexual diaspores on oceanic islands recently pointed to 
a global loss of LDD ability in oceanic bryophyte floras [51]. It 
must be acknowledged, however, that spore production may be 
overlooked when fieldwork has not been sufficiently intensive (e.g. 
[59]). The ability of taxa to undergo LDD via asexual propagules 
(usually larger than 50 |J,m) is still not clearly understood and, 
contrary to classical studies with spores [60], [61], little is known 



about the maximum distances that vegetative propagules might 
actually travel (cf. [62]). Simulation studies have recently pointed 
out that above a diameter of 20 nm wind dispersal of microbes 
between continents becomes increasingly unlikely, and it does not 
occur at all for those of 60 |j,m diameter [63]. Hence, once 
airborne, bryophyte spores can be virtually transported large 
distances by wind, so other establishment impediments such as 
edge colonization or gene surfing [64] and niche specialization 
[65] could be as limiting factors as dispersal per se. 

Here we are assuming that LDD is likely achieved by spores 
because they are resilient, microscopic and released in several 
millions, with some of them eventually being able to colonize an 
island in the very long range. This mechanism might particularly 
apply to archipelagos such as the Macaronesian islands [36], 
which are relatively less isolated and exhibited a higher 
connectivity with the continental sources in the past, due to the 
presence of a higher number of emerged islands that remain today 
as submersed seamounts [66]. Our results are in line with these 
arguments showing that the limited contribution of the Equilib- 
rium Model of Island Biogeography in Macaronesian bryophytes 
can be justified at least by the negligible effect of geographic 
isolation we expected. Different studies have shown similar results 
in bryophytes [17], [67], [68], as well as in other organisms that 
disperse passively by spores, for which wind connectivity seems to 
be more important than geographic proximity between land 
masses [69]. 

By contrast, the effect of area on bryophyte species richness has 
been either supported or not in different studies, both in the case of 
islands and isolated patches on fragmented landscapes (e.g. [22], 
[37], [67], [70]). Our results show that the General Dynamic 
Model -which is ultimately an extension of that originally 
proposed by the EMIB- still exerts significant predictive power 
over island richness after accounting for other factors, particularly 
in the case of mosses (cf. [17]). The variation of island area 
through time is not expected to affect the speciation process in this 
taxon because of its high dispersal potential, but both immigration 
and extinction are probably influenced by the changes in habitat 
diversity through the island's ontogeny, as predicted by this theory 
[8] . The few studies that have evaluated the effects of island age on 
whole floras of spore-dispersed plants show a comparatively lower 
predictive power of this variable than area and habitat diversity 
[17], [29], [67]. At this point, however, one could argue that it is 
very difficult to disentangle the effect of HD and GDM hypotheses 
because time is implicitly accounting for the changes in island 
topography which, in turn, is also correlated with area. In fact, we 
observe strong shared effects between both hypotheses although, 
as we mentioned above, this could be related with the variables 
chosen to represent habitat diversity. Hence, what we may 
interpret from our results is that neither time nor area (nor 
topography) alone explain as much variability in the data as both 
factors together (around 40% or even 60% in liverworts and 
mosses, respectively). 

In spite of its importance, the precise relationship between area 
and species richness is still under debate. Several alternative 
hypotheses have been proposed to explain the importance of this 
relationship on extinction, colonization and speciation rates or 
even stochastic processes [71], among which its correlation with 
habitat diversity seems to apply in different taxa (see [9] and 
references therein). Yet, separating the effects of these factors is 
statistically challenging, and comparatively fewer studies have 
analyzed the influence of the 'effective area' over species richness, 
that is, the relevance of area of suitable habitats [18], [42]. 
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Some notes of caution 

Determining which predictors are the most ecologically 
meaningful for a particular taxon is always difficult to scrutinize 
by the available statistical techniques, particularly because the 
obtained relationships depend basically on (i) the geographic 
extent of the study region, (ii) the dataset (i.e. sample size) and (iii) 
the existence of multicollinearity among predictors. These 
problems are exacerbated in the case of island biogeography due 
to the combination of small sample sizes and different relationships 
between species richness and island characteristics among different 
archipelagos [21]. However, our objective here was to choose the 
best combination of variables to represent each hypothesis - i.e. to 
account for data variation in species richness - rather than 
selecting the individual predictors that are most biologically 
important. In spite of this, we mentioned this issue and, by 
comparing also with analyses including the Cape Verde archipel- 
ago (File SI), we could say that at least orographic mist layer 
together with temperature seem to be important for these 
organisms (see Tables 1 and 2, and Table Sl.l and SI. 2 in File 
SI). Few studies have proved experimentally the influence of mist 
layer because quantitative data on air humidity are often hard to 
obtain (e.g. [46], [72]), hence using typically indirect measures like 
bryophyte cover itself [47]. Further investigation is however 
required to confirm the role of mist precipitation by incorporating 
an actual proxy (e.g. frequency or volume of mist precipitation), 
testing its independent and combined effect with annual precip- 
itation. 

Obviously, widening the spatial extent of the study entails an 
overall higher contribution of climate - in both groups of 
bryophytes - due to a stronger latitudinal gradient (Table Sl.l, 
SI. 2 and Fig. Sl.l in File SI). Nevertheless, habitat variables have 
a notably contribution explaining the richness of mosses (Table 
Sl.l, SI. 2 in File SI), while showing a high shared variance with 
climate (Fig. SI in File SI). Nonetheless, the consistency between 
the results obtained with ordinary least squares models and spatial 
autoregressive regressions indicate that differences in the spatial 
positioning of islands are not conditioning the main patterns found 
at the Macaronesian extent. 

Conclusions 

Our findings indicate that large-scale variations in the 
Macaronesian bryophyte diversity are highly influenced by 
environmental factors but also, at least in the case of mosses, by 
factors related with the island ontogeny. This could be the case for 
other taxa with high dispersal ability. We have also shown strongly 
different macroecological patterns between mosses and liverworts, 
reinforcing the idea that not only dispersal ability, but also 
different ecophysiological responses of these two evolutionarily 
distinct lineages, are probably shaping the distribution of species 
diversity. Our results evidence the seeming importance of climate, 
in particular orographic mist layer, for liverwort diversity, while in 
mosses this factor has a similar or, at least, indiscernible effect to 
that of habitat or even the geologic ontogeny. These results point 
to a presumably large relevance of the establishment process on 
the island diversity of spore-dispersed plants. Future studies, using 
broader spatial extents are required to generalize these conclu- 
sions. 

Methods 

Area of study 

The five Macaronesian archipelagos (Azores, Madeira, Selva- 
gens, Canaries and Cape Verde) lie in the North Adantic Ocean, 



covering a maximum latitudinal extension of almost 3000 km. 
Although all these archipelagos have a volcanic origin, the 
geographical characteristics of the islands vary widely in both size 
(ranging from 3 km 2 to more than 2000 km 2 ) and isolation 
(oscillating from less than 100 km to about 1800 km in distance to 
the nearest continent). Their maximum geological age also differs 
significandy among islands, from the youngest island of Pico in the 
Azores (less than half a million years) to the Selvagem Grande in 
the Selvagens archipelago that dates back to the Miocene (27 Ma). 
It is known, however, that a much older and interconnected 
"Palaeo-Macaronesia" existed during the Paleocene (60 Ma), most 
of which remains today as submersed seamounts [66]. There are 
also evident climatic differences between archipelagos along the 
large latitudinal gradient they form, from the temperate oceanic 
conditions of the Azores to the Mediterranean climate of Madeira, 
Selvagens and Canary Islands. The most extreme conditions for 
bryophyte survival occur in Cape Verde, where tropical arid 
climate prevails, and in the Canarian islands of Lanzarote and 
Fuerteventura, all showing desert affinities due to the Sub-Saharan 
influence. Despite such contrasting island features, there are some 
biotic elements shared between most archipelagos among which 
the evergreen laurel forests (or laurisilva) are probably the best 
representative example for Azores, Madeira and the Canaries. 
Although most of these forest areas were highly reduced after the 
Pleistocene glaciations and current human activity, they present 
the optimal habitat conditions for attaining maximum levels of 
bryophyte diversity (e.g. [73]). 

Data compilation 

We calculated total species richness (S TO t) per island using 
recent checklists updated with some relevant references (see 
further details in [22]). As the species checklist from Selvagens and 
some Cape Verde islands cannot be considered reliable [22], we 
focus our analyses on the main islands (n = 19) of Azores, Madeira 
and the Canaries (but see also File SI, where we show additional 
analyses including some Cape Verde islands that could be 
comparable in terms of inventories). We also run separate analyses 
using the species richness of the two dominant groups of mosses 
(Sm) and liverworts (S£j because they normally present different 
ecophysiological responses [40] . In total, our database included all 
the 729 bryophyte species recorded in the four Macaronesian 
archipelagos. Out of these, 505 are mosses, 218 liverworts and 6 
hornworts. 

Fifteen predictor variables representing geography, time, 
habitat diversity and climate were used to evaluate the relevance 
of the four hypotheses formulated above (see File S2 for further 
details on computation and data sources). For Hi (i.e. EMIB), we 
compiled data on island area (A), distance to mainland (D M ) and 
distance to the closest island [Dj). We also calculated the neighbour 
index (jV) for all the islands as proposed by Kalmar and Currie 
[1 1], to account for the combined effects of the area and distance 
of nearby islands. For H 2 (i.e. GDM), apart from area we obtained 
the time elapsed since island formation (7j for each island. For H 3 
(i.e. HD), we used the number of main ecological zones in the 
islands (EQ as well as three topographical surrogates: maximum 
elevation (ELEV), standard deviation of elevation (sdELEV) and 
diversity of slopes (SLOPEdiv). Finally, for H 4 (i.e. CLIMATE), we 
used six variables accounting for extreme, average and intra- 
annual variation of precipitation and temperature that are of 
particular importance for bryophyte distribution: maximum 
temperature of warmest month (T A1A x), precipitation of driest 
quarter (Pmin), temperature seasonality (Tj), precipitation season- 
ality (Ps), annual precipitation (Pamn) an <f an index of horizontal 
precipitation as surrogate of orographic mist layer (MistE). 
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Statistical analyses 

We used ordinary least squares (OLS) regressions to evaluate 
the influence of the different predictors on bryophyte species 
richness. Prior to the analyses, all these variables were standard- 
ized to zero mean and one standard deviation to remove the effect 
of different measurement scales. We first explored the univariate 
relationships between the considered dependent variables (Stot, 
Sm and Sj) and each one of the 15 predictors. We evaluated both 
linear and curvilinear (quadratic) relationships selecting for 
subsequent analyses the function that maximized the explained 
variance in the response variable. In the case of GDM we included 
both the linear (Tj and quadratic (TT 2 ) functions of time as 
suggested in the literature [8], [10]. We then selected the best 
subset of predictors representing each hypothesis through multiple 
regressions, using the Akaike information criterion corrected for 
small sample sizes (AICc) to compare alternative models. Lower 
AIC C values indicate a compromise between higher model fit and 
lower complexity, so for each hypothesis we chose as best model 
the one with the minimum AICc [74]. In the case of GDM, we 
applied direcdy the ATT 2 model proposed by Whittaker et al. 
(2008) that assumes an unimodal response of species richness to 
time, also accounting for the positive and monotonic influence of 
area. We additionally reported conventional statistics (E 2 and F- 
test) for all the obtained models. 

To ascertain whether the combined effects of the different 
models selected previously contribute to increase significandy the 
explained variance, we made pairwise comparisons between the 
best models chosen for each hypothesis. To separate the single and 
combined effects of the different models, we used partial regression 
analyses [75] . Unlike standard regression methods, this technique 
allows disentangling the proportion of explained variance that can 
be attributed exclusively to one set of factors once the effect of 
other sets has been controlled for, assuming that combined effects 
reflect the shared variance that cannot be unequivocally attributed 
to any of the individual sets of predictors. 

In order to examine whether archipelago' idiosyncrasies non- 
strictiy related with the abiotic island characteristics here 
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